Characterization of dynamical phase transitions in quantum jump trajectories beyond 

the properties of the stationary state 



Igor Lesanovsky, 1 Merlijn van Horssen, 2 Madalin Guta, 2 and Juan P. Garrahan 1 

1 School of Physics and Astronomy, University of Nottingham, Nottingham, NG7 2RD, UK 

2 School of Mathematical Sciences, University of Nottingham, Nottingham, NG7 2RD, UK 

(Dated: November 9, 2012) 

We describe how to characterize dynamical phase transitions in open quantum systems from a 
purely dynamical perspective, namely, through the statistical behavior of quantum jump trajecto- 
ries. This approach goes beyond considering only properties of the steady state. While in small 
quantum systems dynamical transitions can only occur trivially at limiting values of the controlling 
parameters, in many-body systems they arise as collective phenomena and within this perspective 
they are reminiscent of thermodynamic phase transitions. We illustrate this in open models of in- 
creasing complexity: a three-level system, a dissipative version of the quantum Ising model, and the 
micromaser. In these examples dynamical transitions are accompanied by clear changes in static 
behavior. This is however not always the case, and in general dynamical phase behavior needs to be 
uncovered by observables which are strictly dynamical, e.g. dynamical counting fields. We demon- 
strate this via the example of a class of models of dissipative quantum glasses, whose dynamics can 
vary widely despite having identical (and trivial) stationary states. 



Recent experimental progress in quantum optics and 
cold atomic physics has stimulated a great interest in the 
study of open many-body quantum systems [IHZ]- Cur- 
rently, much effort is dedicated to the understanding and 
classification of dynamical phases and transitions among 
them, for example in the Dicke Model [7HS], in lattice 
bosons subject to engineered dissipation pQ, and in spin 
systems [3H51 HH] . Typically, dynamical phase transitions 
are detected and analyzed through changes in static order 
parameters, such as superfluid density or spin polariza- 
tion, calculated in the system's stationary state. 

The aim of this work is to characterize dynami- 
cal phases exclusively through time-correlations within 
quantum jump trajectories, and not through static order 
parameters. We do so by building on an elegant con- 
nection between open quantum systems - described by 
a Lindblad master equation - and matrix product states 
[TT] (MPS) put forward in [12]. We pursue the usual dual 
description of open system dynamics: On the one hand 
any individual realization of the dynamics is stochastic. 
In the case of an open quantum system this is repre- 
sented by a stochastic wave function corresponding to a 
specific sequence of quantum jumps, with the whole en- 
semble of these stochastic trajectories being encoded in 
a MPS. On the other hand, the evolution of probabili- 
ties is deterministic and is derived from the evolution of 
the density matrix under the action of a quantum mas- 
ter operator (QMO). The dynamical phase structure of 
an open system is given by the low lying spectrum of 
the QMO jS]- More specifically, dynamical phases are 
characterized by how the spectrum responds to changes 
in physical parameters [3H5] or counting fields [T31 [H] . 
When this response is non-analytic we have the signature 
of a dynamical phase transition jT3l [14] . This analytic 
structure is, however, also mirrored in the ensemble of 
trajectories. So in order to understand dynamical phase 



transitions we need to study both the spectrum of the 
QMO and the MPS that encodes the trajectories. 

We illustrate our approach through a number of ex- 
amples exhibiting a wide variety of dynamical features. 
To set the stage we first discuss three well-studied open 
systems that would appear to exhaust the range of possi- 
bilities: a "blinking" three-level system [21 [TS], a dis- 
sipative quantum Ising model with a dynamical first- 
order transition [2l 0], and the micromaser, which has 
dynamical transitions of both first and second order kind 
[Trjl [T7] . In all these cases transitions in dynamics are re- 
lated to changes in static behavior. We consider however 
a fourth class of problems, illustrated via a model of dis- 
sipative quantum glasses |10] , whose dynamical behavior 
can change abruptly while their statics remain invariant 
throughout. This highlights the need for a dynamical 
approach like the one presented here for characterizing 
complex open systems where static order parameters are 
non-existent or difficult to identify. 

The density matrix p of the open quantum systems 
we consider evolves under the Master equation dtp = 
W(p) with the QMO given by W(») = %{•) + £>(•). 
Here the super-operators %{•) = — i [H, •] and £>(•) = 

Em=i L ™ • L L ~ \ Em=i { L L L m, •} govern the coher- 
ent and dissipative dynamics, respectively. They depend 
on the hamiltonian H and the set {L m , m = 1, . . . , N} of 
jump operators. The specific form of the latter depends 
on the coupling of the system to the bath. It is well estab- 
lished that an open system dynamics generates a MPS on 
the bath degrees of freedom [T3] . We use this connection 
to link the analysis of the emission sequence of the bath 
quanta, i.e. the quantum jump trajectories, to the more 
familiar picture of static phases of the ground state of a 
one-dimensional spin system. The prescription of Refs. 
[12] is most clear for the evolution of the density matrix 
p over short but finite time intervals 8t, represented for 
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where |^(/)) is the (unnormalized) MPS 



FIG. 1. a: Kraus operators defining open dynamics: Kq is 
the non-unitary (no-jump) evolution, while K m >o refers to 
the occurrence of a quantum jump of kind m. b: The am- 
plitude for a sequence of quantum jumps is obtained directly 
from the Kraus operators. This can be mapped to the state 
of a fictitious one-dimensional spin system, and all possible 
trajectories with their amplitudes can be gathered in an MPS 
interpreted as the quantum state of the spin chain, c: The 
two-time correlation time of quantum jumps (or correlation 
length of the spin chain) is the inverse of the spectral gap 
of the QMO W; dynamical transitions occur when this gap 
closes as a function of an external parameter. 



instance by the Kraus map 



T St (p) = e nst e vst (p) = K Q pKl + V K m pK^, (1) 



m— 1 



with Kraus operators Kq = e 



-iStH 



5tL m . Such discrete (but non- 
unique) representations of dissipative evolutions have 
been employed recently in experiment to simulate the 
dynamics of open many-body systems |18j . 

Fig. [T^, sketches the action of the Kraus operators: 
Kq corresponds to non-unitary (i.e. no-jump) evolution, 
while K m> Q represents the effect of the quantum jump 
associated to L m . The evolution of the density matrix 
over macroscopic times is generated by multiple appli- 
cations of the map ([!]). This produces quantum jump 
trajectories as shown in the upper panel of Fig. [TJd. The 
probability for a certain trajectory {n\,n%, . . . ,um} to 
occur after M time steps is given by p ni n 2 ...n M = 2/ I (/ I 
K nu . . . K n2 K ni | i)\ 2 , where \i) is the initial state of the 
system. The sum runs over a basis of final system states 
| /) , and each term is the probability for connecting the 
initial and final states via a certain sequence of quantum 
jumps. 

These probabilities are encoded in an MPS which can 
be thought of as being generated by letting the system 
interact sequentially with a chain of (JV+ l)-dimensional 
spins initially prepared in a fixed pure state [T2] . After M 
steps the quantum state of the system and the M spins 
with which it has interacted is |\&) = J2f I/) ® IV'Cf))) 



and K„ 



iStH 



N 

W))= E (f\K nM ...K ni \i)\ ni ,...,n M ), (2) 

n M ,...,ni=0 

with the sum running over all spin basis vectors. We 
can think of | ■*/>(/)) as the ground state of a fictitious 
one-dimensional spin system with specific boundary con- 
ditions (see Fig.[l]D). The state |'F) therefore encodes the 
whole ensemble of quantum trajectories: each basis state 
|n.i, Um) corresponds to a specific trajectory and its 
amplitude (/| K„ M . . . K ni \i) is related to the probabil- 
ity Pn 1 n 2 ...n M °f ^ occurring dynamically. While this con- 
nection is only formal it illustrates that the study of dy- 
namical phases of open systems is not different from that 
of the static properties of a one-dimensional spin system, 
regardless of the actual spatial dimension of the open 
problem. Dynamical phase transitions will then become 
visible in the time correlations of the quantum jumps 
which correspond to spatial correlations in the spins. 

The limit of very long times is the "thermodynamic 
limit" of the associated spin problem. In this regime 
the two-time (connected) correlations of observables at 
positions y and y + x have the asymptotic form (see 
e.g. Ref. US]), mf)\A(y ) B(y+^\i;(f)) c cx Re/4 = 
e~ x ^ cos(x<fo)- Here p>2 = \P2\ exp(±i</> 2 ) denotes the 
eigenvalue(s) with the second largest absolute value, of 
the transfer operator E = J2 n =o ® (or equiva- 
lently of Tg t ), the largest eigenvalue being p,\ = 1 by the 
conservation of probability of the map 0. The corre- 
lation length £ is given by = — log | /j,2 | ■ The cor- 
relations exhibit exponentially damped oscillations when 
the eigenvalue P2 is complex, in which case p* 2 is also 
an eigenvalue of T$ t - It is more convenient to directly 
study the eigenvalues A of the QMO W. In the limit 
St — > temporal correlations between quantum jumps 
behave as (A(t)B(t + t')) c cx exp(-i'/r) cos (uit'), with 
t = — l/Re(A 2 ) and ui = Im(A 2 ), where A 2 is the eigen- 
value of W with the second largest real part, Fig. [TJ; 
(note, the eigenvalue with largest real part is Ai =0). 

Dynamical phase transitions will manifest in the clos- 
ing of the spectral gap of the operator W, i.e. A 2 — > 
(see Fig. [TJ;), a well-known feature in statistical mechan- 
ics which has also recently been used to characterize dy- 
namical transitions in quantum spin models [S] and the 
micromaser [17] . This closing of the gap can occur in 
qualitatively different ways giving rise to a different kinds 
of dynamical transitions. In the following we discuss this 
with four examples. 

(i) Three-level system - The three level system [15] 
as depicted in Fig. [2^, is described by the Hamiltonian 
H 3 = Qi |0)(1| -I- f2 2 |0)(2| + h.c, and a single jump op- 
erator L = v / k|0)(1|. This system exhibits a dynamical 
phase transition at f2 2 = (where A 2 =0), which mani- 
fests itself in a strongly intermittent behavior of photon 
emission when f2 2 <C [2D]. The physical reason for this 



3 




FIG. 2. Sketch of model systems, real part of the spectrum of the QMO, and sample trajectories (from top to bottom), a: The 
driven three-level system has a dynamical transition at the limiting value Q.2 = (circled region). For £^2 > quantum jumps 
are highly intermittent (we show the case ^2 = 0.01k and Q.\ — 4k). In the boxed region A2 acquires an imaginary part and 
time-correlations become oscillatory, b: The open Ising model with L spins has a genuine many-body dynamical first-order 
transition and the spectral gap closes only in when L — > 00. In the parameter regime chosen (V = 100k) the spectral gap 
becomes small for 10k < Q, < 30k. In this near coexistence region of two dynamical phases one observes strongly intermittent 
behavior of quantum jumps (we show data for L = 10, f2 = 25k). c: The micromaser displays both first and second-order 
dynamical transitions. The two sample trajectories are taken at the second-order transition point (red circle/trajectory) and 
within the coexistence region of two dynamical phases (blue box/trajectory). In order to facilitate the representation we have 
collected quantum jumps in time bins of length ibin = 3/k. 



is that at f^2 = the system decouples into a driven two- 
level system undergoing frequent quantum jumps and an 
inactive single dark level. This corresponds to a twofold 
degeneracy of the leading eigenvalue (Fig. [2£i) , which is 
lifted when > 0. For ^ the system can switch 
between these two phases on a timescale fi^" , which re- 
sults in strongly intermittent behavior reminiscent of a 
(smoothed) first-order transition |14j . From the perspec- 
tive of the MPS this corresponds to the quantum phase 
transitions reported in Ref. [12]. Beyond the transition 
at f^2 = 0, the three-level system features also a dynami- 
cal transition at finite ^2 where time correlations become 
oscillatory (see Fig. [2^,). Related transitions have been 
reported in an NMR experiment; see Ref. pi] [22j. 

(ii) Dissipative Ising model - In a many-body system 
the degeneracy leading to a phase transition does not 
need to be imposed externally, but can appear sponta- 
neously in the thermodynamic limit. An example is the 
dissipative quantum Ising model [SIS] sketched in Fig.J^. 
The Hamiltonian is that of a one-dimensional Ising chain 
in a transverse field: H \ = il J2k=i °x + ^ J2k=i a z a z +1 
where a\ are the Pauli spin matrices. The jump opera- 
tors are = \[k.o^ , which produce incoherent spin flips, 
It) — ^ \i) > at a rate n. In Fig. ^jp we show the real part of 
the eigenvalue spectrum of the QMO W for various sys- 
tem sizes L. With increasing L the spectral gap closes 
over an entire (coexistence) region in parameter space 
and we expect it to approach zero when L — > 00 and 
hence the dimension of the transfer operator E becomes 
infinite. This resembles a phase transition in the ther- 
modynamical sense within the fictitious spin system. For 



finite L and finite gap the system switches between two 
dynamical phases on long but finite timescales and quan- 
tum jump trajectories are strongly intermittent. Like in 
the three-level system, in this case the dynamical tran- 
sition can be traced back to a bistable static behavior: 
When the emission of photons is plentiful the average 
magnetization is close to zero, while it is large and neg- 
ative during the dark periods [4]. 

(Hi) Micromaser - The micromaser features a critical 
point and a sequence of first order transitions [El [17] . 
It is modeled by a resonant single-mode cavity coupled 
to a finite temperature bath and pumped by excited 2- 
level atoms which are sent into the cavity with a constant 
rate r, Fig. [2j:. The Hamiltonian is zero and there are 
four jump operators, two from the atom-cavity interac- 
tion, L\ = \[ra^ sln (^jj_ a ) anc [ £ 2 _ y/r cos ^(f)V aa^ , 

and two from the cavity-bath interaction, L3 = t/ko. and 
L4 = ^Jva) . Here a, are the raising/lowering operators 
of the cavity mode, k and v are the thermal relaxation 
and excitation rates, and (j> encodes the atom-cavity in- 
teraction [T7]. The spectrum of W is real and is de- 
picted in Fig. [2J: as a function of the "pump" parameter 
a — (f>/\/r. A sequence of first order transitions are vis- 
ible beyond a — 4 and quantum jump trajectories (here 
we monitor quantum jumps associated to L{) show the 
typical intermittent behavior. In contrast, in the vicinity 
of a = 1 the spectral gap closes in a way that makes the 
spectrum dense. This is the onset of a second-order phase 
transition which strictly only occurs in the limit r — > 00 
|23) . Typical quantum jump trajectories near the critical 
point fluctuate very strongly, as shown in Fig. [2j;. Also 
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FIG. 3. Dissipative quantum glass model, with glassiness con- 
trolled by the parameter p: at p = system is non-interacting, 
for p — 1 dynamics is fully constrained. The spectral gap 
closes at p = 1 (we show the case of k = 8f2). Trajecto- 
ries shown consist of 5000 quantum jumps each. Dynamics 
changes drastically with p but the stationary state remains 
invariant. Close to p = 1 emission periods are localized in 
time and space (see magnified region). For ease of visibility 
jumps are collected in 500 evenly spaced time bins. 



the micromaser dynamical transitions are accompanied 
by a change in the statics, i.e. at the transition points 
the photon occupation of the cavity undergoes a sudden 
change: at first-order transitions the mean photon num- 
ber switches between two largely distinct values, while 
at the critical point the variance of the photon number 
undergoes a jump [IB]. This seems to suggest that dy- 
namical transitions can always be anticipated to occur 
by simply considering static or steady state properties. 
The next example proves that this is not the case. 

(iv) Dissipative quantum glass - This system is related 
to the dissipative quantum glass models of [10]. It is a 
spin chain with Hamiltonian H g = ^Efci ^i/t+iW 
and jump operators Lk = y/K&k fk+i(j>) (see Fig. 
[3]). The operators fhip) a- re kinetic constraints. If 
fk(p) = 1 the system is just a set of non-interacting 
two level systems with steady-state density matrix p ss = 



®k=i { 2 + k 2 +Lj 2 } p k + { 2 re 2+L 2 } 



where ui 



Here the statics is clearly featureless. The problem be- 
comes interacting and glassy if we choose the constraints 
to be f k {p) = pQ k + (1 -p)l with < p < 1: For p = 
we have the non-interacting problem and for p = 1 a 
fully constrained quantum glass |10j . When p = 1 the 
state of the spin on site k can only change if the state 
\(p) k+ i of its neighbor satisfies Qk+i \<P)k+i ^ 0, leading 
to correlated dynamics in the system. The parameter p 
controls how glassy this dynamics is. Fig. [3] shows that 
the spectral gap of the QMO would close at p = 1. Here 
we expect a dynamical first-order phase transition to oc- 
cur in analogy with classical constrained models [131 [21] • 
A crucial feature of this model is that the steady state 
for any value of p is the trivial p ss , i.e., there is no change 
in the statics despite the change in the dynamics, as is 
evident from the example jump trajectories shown Fig. 
[3] As p — >• 1 the system is most of the time inactive 
and quantum jumps become more and more localized in 



space and time, a phenomenon called dynamical hetero- 
geneity which is a hallmark of glassy relaxation [25] • In 
contrast to the three systems described above, here the 
dramatic dynamical change with p is impossible to guess 
from static properties which are p independent and trivial 
throughout. 

Since the ensemble of trajectories is fully encoded in 
the MPS 1^) the sample trajectories of Fig. [3] indicate a 
transition within the MPS as p — > 1. Usually such static 
quantum transitions are accompanied by a singularity, 
such as a logarithmic divergence in the entanglement en- 
tropy of large spin blocks [3BJ . In our case this would be 
an entanglement between quantum jumps in subsequent 
long time segments. The entropy of a large block, how- 
ever, is just two times the von Neumann entropy of the 
stationary state Se = — 2Tr(p ss log p ss ) of the system, 
where p ss = Tr^|5')(^I'| |27j . Hence, due to the invari- 
ance of the stationary state with p, this entanglement 
measure does not detect the transition. 

The reason for this is that in this problem — as is 
generic in both classical and quantum glasses, and likely 
to occur often in complex many-body systems — the ap- 
propriate fields driving the transition do not couple in an 
obvious way to static quantities, but do couple directly to 
time-integrated observables. An instance of these are the 
"counting" fields introduced when computing full count- 
ing statistics [28] of dynamical observables [T3J [14] . Con- 
strained models such as the one described here are known 
to exhibit dynamical transitions in trajectories, which 
arc evident in the moment generating function (MGF) of 
the number of quantum jumps, Z t {s) = J2jPt(J)e sJ , 
where Pt{J) is the probability of observing J jumps in 
time t. In the t — > 00 limit Z t (s) becomes singular at 
some value s = s c of the counting field [131 HI], indicating 
a phase transition in the ensemble of trajectories. At long 
times, the MGF is obtained from the largest eigenvalue 
of a deformation W —> W s of the QMO parametrized 
by s [Ml HI]. The field s therefore couples directly to 
the spectrum, so that a singularity of the MGF at s c in- 
dicates the existence of close to degenerate but distinct 
dynamical states. By driving s one can single these states 
out; see Refs. gHHUT] for details. 

The MGF can be connected to the MPS through the 
norm, Z t (s) = (*(s)|*(s)), of the deformation |*(s)) = 
e -sJ/2| v j / ^ wriere j j s an operator that counts the num- 
ber of n m> o in a state |m, Hm). In a thermodynamic 
analogy Z t (s) is like a partition sum over trajectories and 
s a chemical potential which favors/disfavors quantum 
jumps. The state |^(s)) is thus a superposition of MPS 
like those of ([2| but with each term weighed by a factor 
of e~ s / 2 for each jump. The eigenstate p(s) of W s cor- 
responding to the largest eigenvalue is related to |^(s)) 
through T)ty|*00)(*O)| = Z t (s)p(s), where Trp(s) = 1 
and p(0) = p ss . The entanglement entropy of this state, 
Se = — 2Tr[p(s) logp(s)], will depend on s, and will dis- 
play non-analytic behavior at s c . At this level counting 
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fields work in a similar manner to more standard static 
fields which drive phase transitions, but couple directly 
to the relevant dynamical order parameters that reveal 
transitions in quantum jump trajectories. This perspec- 
tive should be useful in the study of dynamical phase 
transitions in systems where they are not obviously con- 
nected to a change in spatial correlations. 
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